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Abstract. A simple proof of a weak limit theorem for transition probabilities of quantum walks 
on the one-dimensional integer lattice is given. To prove it, an abstract algebraic structure for 
one-dimensional quantum walks is introduced. This structure characterizes, in some sense, one- 
dimensional quantum walks. A natural computation using this algebraic structure leads us to obtain 
an effective formula for the characteristic function of the transition probability. Then, asymptotic 
behavior is deduced by using simple properties of the Chebyshev polynomials. 
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1. Introduction 



< 

■ The notion of quantum walks was introduced by Aharonov-Davidovich-Zagury (PQ) as a quan- 

tum analogue of classical random walks and then re-discovered in computer science. In particular, 
Ambainis-Kempe-Rivosh ([2]) utilized two-dimensional quantum walks to improve Grover's quan- 
tum search algorithm. See [5], [7] for historical backgrounds and various aspects of quantum walks. 
Among numerous works, Konno ([6]) obtained a weak limit formula for transition probabilities of 
quantum walks on the one-dimensional integer lattice. In this short paper, we present a simple 
proof of Konno's weak limit theorem. First of all, let us state Konno's theorem. The quantum 
walks we consider in this paper is defined in terms of a two-by-two (special) unitary matrix 

§: (i) A= {-1 a)> a ' 6GC > H 2 + I6| a = i. 

, and its decomposition, 

Let £ 2 (Z) <8> C 2 be the Hilbert space of square summable functions on Z with values in C 2 whose 
inner product is given by 

{f,9)=J2( f( x )i9{x) )& (/, 9 G ^) ® c 2 ), 

where (-,-)c 2 denotes the standard inner product on C 2 . For any u G C 2 and x E Z, define 

5 x ®uel 2 (Z)®C 2 by 

u (y = x), 



(5 X (8) u)(y) 



(y^x). 



Then, the vectors, 5 X ® e, {i = 1, 2, x G Z), where {ei,e2} denotes the standard orthonormal 
basis in C 2 , form a complete orthonormal system in £ 2 (Z) (%> C 2 . The unitary evolution, U, for the 
quantum walk is a unitary operator on £ 2 (Z) ® C 2 defined as 

(2) U = Pt + Qt-\ 



Date: October 3, 2012. 

The author is partially supported by JSPS Grant-in- Aid for Scientific Research (No. 21740117). 

1 



2 



TATSUYA TATE 



where r is the shift operator on £ 2 (Z) ® C 2 defined by t(5 x ® u) = 5 x +i ® u. Since the operator U 
is unitary, the function 



defines a probability distribution on Z for each positive integer n and a unit vector tp in C 2 . We 
call the distribution p n ((p;x) the transition probability of the quantum walk In this context, 
Konno's weak limit theorem is stated as follows. 

Theorem 1. Let ip = t (p>i,<P2) be a unit vector in C 2 . Suppose that the components, a, b, of the 
matrix A given in (pTJ) are non-zero. Then, we have the following weak limit formula. 



There are several methods to prove Theorem [TJ In [6], Konno computed the probability distri- 
bution p n (p\x) directly and explicitly, and then used a result on asymptotic behavior of Jacobi 
polynomials ([3])- In [3], Grimmett-Janson-Scudo employs an integral formula involving the eigen- 
values of the matrix 



Explicit computation of p n ((p;x) in [6] involves computation of words in the matrix P, Q. It is a 
very interesting computation, but it is a bit complicated. Furthermore, the asymptotic formula for 
Jacobi polynomials used in [5] is not elementary. The integral formula involving the eigenvalues of 
the matrix A(z) is very powerful. Indeed, in [9], the integral formula is used to compute the local 
asymptotics of p n (p;x). However, the integral formula does not make explicit form of p n (ip;x) 
completely clear. Furthermore, it would not be quite easy to compute the eigenvalues in higher 
dimensional cases. 

In the present paper, we give a simple and systematic computation of p n (ip; x) and prove Theorem 
[I] in an elementary way. Our strategy is to exploit an algebraic structure behind the quantum walk 
U. See Section [2] for details on this algebraic structure. This structure characterizes, in some sense, 
the quantum walk U. Furthermore, this structure naturally leads us to obtain an effective formula 
of p n ((p;x) in a systematic way, and then asymptotic behavior of the characteristic function of the 
distribution p n (ip;x) is easily deduced using simple properties of the Chebyshev polynomials. In 
the present paper, higher dimensional cases are not discussed. (See [10] for a pioneer work on the 
weak limit formula of certain one-parameter family of two-dimensional quantum walks.) However, 
we expect that the algebraic structure introduced here would give some insight for quantum walks 
in higher dimension. 



*In [6], the 'row' decomposition of the given matrix A is used. In [9], the unitary operator Pt + Qt is used. In 
each case, suitable change of the initial data or the change of the sign of the variable x will give us our distribution. 



(3) 






z G C\{0}. 
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2. AN ALGEBRAIC STRUCTURE OF QUANTUM WALKS 
Suppose that we are given unitary operators, 

V, W, a, 

on a Hilbert space T~Lq, whose inner product is denoted by (•,•), such that they satisfy the following 
conditions. 

(QW1) W 2 = -I. 

(QW2) VW = WV~ 1 . 

(QW3) aW + Wo = aV - Va = 0. 

(QW4) a* = a. 

Note that (QW3) implies that a ^ I. Since W, o are unitary operators, (QW1), (QW4) imply 
W* = W' 1 = -W, a* = a- 1 = a. Let 

tt± :=\{I±a) 

be the projection onto the eigenspace of a with the eigenvalue ±1. Since a ^ I, we have 7r_ ^ 
and (QW3) implies that Vir± = ir±V, Wtc± = tt t W. We set 

(4) X = ^(V + V*), Y = -(V-V*), T = X + iaY. 

Then, it is obvious that we have 

XY = YX, XW = WX, YW + WY = 0, VT = TV, TW = WT, 
^ Xa = aX, Ya = aY, Ta = oT. 

Since T* = X — iaY, we see 

rp*rp rprji* X^ ~\~ Y^ I 

which indicates that T is a unitary operator on T~Lq. The following lemma can easily be proved. 

Lemma 2. (1) We have T = -k+V + tt-V*, V = vr+T + vr_T*. 
(2) Set e = VW. Then we have 

e* = e~ 1 = -e, en± = n T e, eW = -V, We = -V* , eV = V*e, ea + ae = 0. 
Example 1. We take ct,/3 G C with \a\ = |/3| = 1 and define the matrices, Vq, Wo, a, by 



a) 1 u \-/3 oy' V° 

The matrices, Vq, Wq, a, are unitary matrices and satisfy the conditions (QW1)-(QW4) with 
T = al, where T is defined in ([!]). 

Example 2. Let V, W be the quantum walks defined by the formula ([2]) with the unitary matrices, 
Vo, Wo, in Example [TJ Then, the unitary operators, V, W, a, on £ 2 (Z) (8) C 2 also satisfy the 
conditions (QW1)-(QW4) with T = ar, where r is, as in Section dj the shift operator. 

Next, we impose the following condition. 

(QWR) There exists a unit vector e £ n+Ho such that ( V x e, e ) = for any positive integer x. 

Note that the conditions (QW1)-(QW4) have the meaning for unitary elements, V, W, a, in a 
C*-algebra, say ^4o- I n this context, the condition (QWR) is a condition on a representation on Ho 
of the C*-subalgebra in Aq generated by V, W, a. 
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Example 3. The unitary operator V on £ 2 (Z) (g> C 2 given in Example [2] satisfies the condition 
(QWR) with the unit vector e = 5q <8> e\. However, the unitary matrix Vo on C 2 given in Example 
[I] does not satisfy the condition (QWR). 

In what follows, we suppose that the operator V satisfies the condition (QWR) with a unit vector 
e G tt+Wo. 

Lemma 3. The vectors T x e (x G Z) form an orthonormal system in fto- 

Proof. It is enough to show that (T x e,e) = for any non-zero integer x. Since V and a are 
commutative, we have 

(6) Te = Ve G ir + n , T~ l e = T*e = V*e G 7r+ft . 

Since T and V are commutative, we see T x e = V x e for any integer x, which combined with the 
assumption (QWR) shows the assertion. □ 

We set 

ef = T x e, ef = T x ee (x G Z), e x = e?, e 2 = ejj. 

Then, by Lemma [2] and the fact that T preserves the eigenspaces of a, we have ef G tt + / Hq, 
e| G TT-Tio, and the vectors (i = 1, 2, x G Z) form an orthonormal system in %q. Now, let £ be 
the closure of the subspace spanned by {ef ; x G Z}. Then, the closure of the subspace spanned by 
{ef ; x G Z} is Since i C ir+Tio, ^ C n-Ho, these are orthogonal to each other. We define the 
subspace 

in fto, in which the set {ef ; i = 1, 2, x G Z} becomes a complete orthonormal system. 

Lemma 4. The operators, V , W , a, and their inverses preserve the subspace ft so that these define 
unitary operators on ft . 

Lemma H] follows from (J6j) and Lemma [2J Indeed, we have the following. 

(7) Vef = e x l +1 , Ve x 2 =e x ~ l , We x = e x+1 , We x 2 = -e x ~ l . 
Our main interest is in the operator, U, defined by 

(8) U = sV + tW, s,t€R, s 2 + t 2 = l. 

It is obvious that U is a unitary operator on fto and on ft. Since the operators of the form ±V, 
±W, a also satisfy the conditions (QW1)-(QW4), (QWR), we may assume that < s,t < 1. Let 
ij) = t (ipi,ip2) be a um t vector in C 2 , and define \f = ipie± + V^^, which is a unit vector in ft. For 
any integer n, we define a function q n (ip; x) in x G Z by 

(9) g n (^x):=KC/^,e?)| 2 + |(C/^,ef)| 2 . 

Since U is a unitary operator on ft, we have Ylx&z 9n(ip] x) = 1 and hence it defines a probability 
distribution on Z. 

Example 4. Let V, W, a be the unitary operators on £ 2 (Z) ® C defined in Example El Let 
e = <5o <8> 6i, which satisfies the condition (QWR) as in Example [3j In this example, the unitary 
operator U defined in ([8]) has the form 
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which coincides with the quantum walk defined in ([2]). Since VW = ( -5 I , we see &2 



-a/3 

yi/F(<5o <g> ei) = — a/3(5o <8> e 2 ). Since T = ar as in ExampleEl the vectors, ef , are given by 

ef = a x (8 x ® ei), ef = -c^ -1 /^ <g> e 2 ) (x G Z). 

Thus, we have "Ho = ~H = £ 2 (Z) <g) C 2 . As above, let V — * (^1 , ^2) be a unit vector in C 2 . Then 
the corresponding vector \I/ in £ 2 (Z) <g> C 2 is 

^ = V'lei + V>2e2 = 5 <S> <p, p = t (ipi, -a/3ip2)- 

From this we have 

(U n *,e x 1 ) = a x (U n (5 ®ip),5 x ®e 1 ), 
( IT 1 *, ef ) = -a~ x+1 l3( U n (5 ®tp),5 x ®e 2 ), 
which implies q n (ip;x) =p n ((p;x) with (tpi,tp 2 ) = {<Pi } — ®f3if2), where p n (tp; x) is defined in ([3]). 
As in Example El it is enough to prove the following theorem to show Theorem (TJ 

Theorem 5. Let tp = *(?/>i, ^2) be a unit vector in C 2 . Let s, t be real numbers satisfying < s, t < 
1, s 2 + t 2 = 1. Then we have 

Wdim^fe^^j^/n = X(- S>S )(y)— 2 ^ T= 2 2 d y> 

n ^°° x& tt(i - r)v s - v 

where the function X(ip; s,t) is given by 

\(i>;s,t) = |V'i| 2 -|^2| 2 + 2Re(Vi^ 2 )-- 

s 

We give a simple proof of Theorem [5] in the following two sections. 

3. Computation of the distribution q n (tp;x) 

The unitary operator, U, defined in ([8]) is written as U = x + iy + w with x = sX, y = sY, 
w = tW, where X, Y are self-adjoint operators defined in @. By ([5|), we see that the operators x 
and iy + w are commutative and 

(10) (iy + w) 2 = ~{y 2 + t 2 ), x 2 + y 2 + t 2 = I. 

Therefore, we get the following. 



U n = (x+ (iy + w)) n = (?) x n ~ k (iy + 

fc=o ^ 7 

[n/2] [(n-l)/2] 

= £( 2 > n - 2 v-/y+ e 

1=0 V 7 1=0 



2l n +1 jx n - 2l -\x 2 -Ly( iy + w ). 

Next, we recall a definition of the Chebyshev polynomials (see [8] for details). The Chebyshev 
polynomials of the first kind, T n (x) (of degree n), and the second kind, f7 n _i(x) (of degree n — 1), 
are defined as 

T n (cos#) = cosnO, U n -i(cos9) = — . 

sin 6 

These polynomials can be written in the following form. 

[n/2] . - [(n-l)/2] 



r„W = £( 2 ")*"-V-i)', E /„- 1 (x)= £ ( 2i " +1 ) 

1=0 V 7 1=0 v 7 
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From this we can write U n as 

(11) U n = T n {x) + U n ^(x)(iy + w) = T n (x) + (iy + w)U n - X {x), 

where x, y and w are operators defined at the beginning of this section. A direct computation 
shows that 

xe\ = |(T + T _1 )ei, xe 2 = |(T + T _1 )e 2 , 
= -(T- r _1 )ei, %e 2 = -|(T- r _1 )e 2; wei = *Te 2 , u;e 2 = -iT _1 ei. 
Let Pn(z), q % n (z) (i = 1,2) be Laurent polynomials (with real coefficients) defined by 
pl(z) = T n {s{z + z- r )/2) + ^(z- z-^Un-Mz + z" 1 )^), 

p 2 n (z)=tzU n . 1 ( S (z + Z- 1 )/2), 

qiiz) = -tz' l U n ^(s(z + z- l )/2), 

q 2 n (z) = T n (s(z + z~ l )/2) - S -{z- z' l )U n ^{s{z + z" 1 )^). 

Since T commutes with all the operators, we have 

U n e 1 =p 1 n (T)e 1 +pl(T)e< 2 , 

U n e 2 = q 1 n (T)e 1 +q 2 n (T)e 2 . 

For any Laurent polynomial p(z) in z £ C, let us denote the coefficient of z x (x 6 Z) in p(z) by 
C x (p). Then, for any x£Z, (jl3j) shows the following. 

(E^e?) =V'iCx(pi)+^2C x (gi), 
(^,61)=^^)+^^). 

Since p l n (z), q l n {z) have real coefficients, we have obtained the following lemma. 
Lemma 6. The probability distribution q n (tjj;x) on Z defined in ([9]) is written as 

q n (ip;x) = \^\ 2 [c x {p l n ) 2 + C x {p 2 n ) 2 ] 

(15) +\H 2 [c x (q 1 n ) 2 + C x (q 2 n ) 2 ] 

+ 2Re(Vi? 2 ) [c x {pl)c x {ql) + C x {p 2 n )c x {q 2 n )} . 
In particular, the characteristic function 

is given by 

(16) E n (0 = |^| 2 [pi(0 + P n 2 (0] + |^2| 2 [Q^(0 + Ql(0] + 2Re (Vi? 2 ) + , 
where the functions -P^(£)> Qn(£) anc ^ ^n(£) = are gu>en fry 

(17) P*(£) = £c^)V^, Q^H^C^) 2 ^, = £ ^(pfjc^e**. 
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4. Proof of the weak limit theorem 
Let p(z), q{z) be two Laurent polynomials, which are, in our notation, written as 

P(z) = J2c x (p)z x , q{z) = ^C x {q)z x . 
Then, we have the following convolution identity. 

(18) y2c x (p)c x (q)w x = / p(wz)q(z- 1 )——. 

We use the formula (fl8|) to deduce the asymptotic properties of P*(£/n), Q l n {^/ n ) an d R l n {£,/n) as 
n tends to infinity. To compute it, we need the following lemmas. 

Lemma 7. For 9 € R, we set f(9) = Cos" 1 (s cos 6>), ui/iere < s < 1. VFe fix ^ £ R. Then, as 
n — > oo, we have the following. 

T n {s cos(9 + = T n (s cos 0) cos(£/'(0)) - sin(n/(0)) sin(£/'(0)) + 0(l/n), 

t7 w -i(acoB(g + e/n)) = ^(acosg) cos(£/'(fl)) + T n (s cos 6) Bm ^ f '^ + 0(l/n), 

sin/(0) 

where 0(l/n) is uniform in 9. 

Proof. Noting that T n (scos(6» + £/n)) = cos(n/(0 + f/n)) and < Cos _1 (s) < f{9) < Cos _1 (-s) < 
7r, a simple application of the Taylor expansion shows the lemma. □ 

Lemma 8. Let k £ Z and £ e R. We se£ 

A,,*(0 = / **T n ( a (e*/ n * + e-^z- 1 )/2)T n ( S (z + z~ 1 )/2) 

J|z|=l 



J3n,fc(0 = / z k T n {s{e i ^z + e-^ n z- 1 )/2)U n ^{s^ + ^ 1 )/^) 

J |^| = l Z7TIZ 

C n , k (0 = [ z k U n ^(s(e^ /n z + e-^z- 1 )/2)T n ( S (z + z" 1 )/2) -^L 

J|z| = l 



£>„,*(£) = / z k U n ^(s(e^ n z + e-^- 1 )/^^^^^ + z~ 1 )/2) 



\z\=i 2niz 
Then, we have the following. 

l + (-l) fe f s , . ( £Vs 2 -x 2 \ dx 



l + (_l)fc rs 
lim A nk (£) = / e fc (x)cos 

n->oo 47T J_ s 

lim B n fe (£) = - lim C n fe (£) 



S 2 — X 2 



n— ^oo 



i-(-i) fe r , , . / evs 2 - x 2 \ dx 

ek{x) sm 



4vr 7_ s I VT^ 1 / ^/(l-x 2 )^ 2 



.r 



21 



1 + /■« 

lim D nk (£) = / e fc (x)cos 

n-voo 47T J_ s 

where we set ek(x) = exp[ikCos~ 1 (x/s)]. 



l + (-l) fc f s , . / £Vs 2 - x 2 \ dx 



Vl — x 2 ) (1 — x 2 )V s 2 — X 2 
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Proof. Since T n {-x) = (-l) n T n (x), C/„_i(-x) = (-l^Un-^x), we see 

\k /*7T 



A n ,k(0 = 1+ i ^ I* e lk9 T n (s cos(6 + t /n))T n (s cos 8) d0, 

27T Jo 

Bn,k(0 = 1 ~ r e lk6 T n (s cos(6 + an))U n ^(s cos 6) d9, 

27r Jo 

Cn,fc(0 = ^i"^ I* e ike U n - 1 {scos(e + i/n))T n (scose)de, 

A^O = r e ike U n ^(s cos(9 + C/n))U n ^(s cos 9) (16. 

27r Jo 



Note that we have 

7T 



/■7T /*7T 

/ |T n (scos0)| d6 < tt, / |J7 n _i(scos6»)| d6 < 
Jo Jo 



vT^ 2 ' 

Thus, by Lemma [7] and the change of the variable, ip = f(0), where f{9) is the function introduced 
in Lemma [71 we have 

l + (-l) fc 



2tt 



CoS " ( ~ s) ( \ 2 / \ ^V^-cos 2 ^ siny 

e fc (cosy?)cos (7199) cos : = 

Cos- 1 ^) V sm<£ / ,/ s 2 _ cos 2 - 



CoS ^"^ / ^s 2 _ cos 2 p\ s i nifi 

efc(cos ip) sm(nip) cos(rnp) sin — ^rfy> + Q(l/n) 

Cos- 1 ^) \ sm^ J ^Js 2 - cos 2 ip 

Since , 1 is an L 1 -function on the integration interval, the Riemann-Lebesgue lemma gives 

Y S 2 — COS 2 If) 



V A m l + (~l) fc [ Cos ~ { ~ s) , , ^s 2 -cos 2 ^ \ siny 
hm Ai,fc(£) = : / e fc (cos^)cos : = dtp 



4vr Jcos-^s) \ sini^ y ^2 _ C os 2 y? 

l + (-l) fc f .. ( ^s 2 -x 2 \ dx 
ej.(x) cos 1 



4vr J_ s \ y/i — x 2 I \/ s 2 — x 2 

Formulas for B n ^(£,), C nj fc(£) and D n ^{C) can be obtained in a similar way. □ 
Now, it is rather easy to prove Theorem [SJ 

Proof of Theorem A direct computation using (|18p and Lemma [8] leads us to the following. 



Pfa/n) = A nfl (0 - S - (fl„,i(0 - Bn,-i(0) + | (e* /n Cn,i(0 - e-«/ B C n ,_i(0 

„2 



v e i5/n ^n, 2 (e) - 2cos(f/n)A M) (0 + e -*/ n I> n) _ 2 (£) 

1 r / CVs 2 -x 2 \ i f ( £Vs 2 -x 2 \ dx 

cos — ; = H — / sin 



2^ 7_ s \ vT^ 2 / Vs 2 - x 2 vr 7„ s I vT^ 2 7 vT^ 2 



1 r / eVs 2 - x 2 \ Vs 2 - x 2 , 
H / cos — ; tt- dx in — > 00). 



Note that the integrand are all even function. Changing the variable y = y \z§ri we have 

hm pi(e/n) = - / {(1 + y 2 )cos(iy) + 2%sin(£ 2 /)}- ^ ■ 

rwoo vr J (1 _ y 2 )^/s 2 - y 2 
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Similar computation shows that 

t f S „s dy 



lim Qli^/n) = lim P^/n) = - f cos(£y) 

n— >oo n— >oo 7T Jq 



s 2 — y 2 



dy 



lim QlWn) = - [ {{l + y 2 ) cos(^y) - 2iy sin(£y)} — 

n^oo 7T J (1 _ y 2 )^s 2 - y 2 

lim «(£/n) + i£(£/n)) = — / y sin(^)- % 

Noting that integrals of odd functions on the interval (— s, s) vanish, we have, by Lemma [6l 

lim E (£/n) - - f c* v l + X ^\ s ^ fUl 
^ En{W - vr J_ s e (l-y 2 )^^ dy ' 

which completes the proof of Theorem [5j □ 
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